Search for: All records

We study the statistical properties of O vi, C iv, and Ne viii absorbers at low-z (i.e. z < 0.5) using Sherwood simulations with ‘WIND’ only and ‘WIND + AGN’ feedback and massive black-II simulation that incorporates both ‘WIND’ i.e. outflows driven by stellar feedback and active galactic nucleus (AGN) feedbacks. For each simulation, by considering a wide range of metagalactic ionizing ultraviolet background (UVB), we show the statistical properties such as distribution functions of column density (N), b-parameter and velocity spread (ΔV90), the relationship between N and b-parameter, and the fraction of Ly α absorbers showing detectable metal lines as a function of N(H i) are influenced by the UVB used. This is because UVB changes the range in density, temperature, and metallicity of gas contributing to a given absorption line. For simulations considered here, we show the difference in some of the predicted distributions between different simulations is similar to the one obtained by varying the UVB for a given simulation. Most of the observed properties of O vi absorbers are roughly matched by Sherwood simulation with ‘WIND + AGN’ feedback when using the UVB with a lower O vi ionization rate. However, this simulation fails to produce observed distributions of C iv and fraction of H i absorbers with detectable metals. Therefore, in order to constrain different feedback processes and/or UVBs, using observed properties of H i and metal ions, it is important to perform simultaneous analysis of various observable parameters.

 

We measure the mean free path ($\lambda _{\rm mfp,H\, \small {I}}$), photoionization rate ($\langle \Gamma _{\rm H\, \small {I}} \rangle$), and neutral fraction ($\langle f_{\rm H\, \small {I}} \rangle$) of hydrogen in 12 redshift bins at 4.85 < z < 6.05 from a large sample of moderate resolution XShooter and ESI QSO absorption spectra. The fluctuations in ionizing radiation field are modelled by post-processing simulations from the Sherwood suite using our new code ‘EXtended reionization based on the Code for Ionization and Temperature Evolution’ (ex-cite). ex-cite uses efficient Octree summation for computing intergalactic medium attenuation and can generate large number of high resolution $\Gamma _{\rm H\, \small {I}}$ fluctuation models. Our simulation with ex-cite shows remarkable agreement with simulations performed with the radiative transfer code Aton and can recover the simulated parameters within 1σ uncertainty. We measure the three parameters by forward-modelling the  Lyα forest and comparing the effective optical depth ($\tau _{\rm eff, H\, \small {I}}$) distribution in simulations and observations. The final uncertainties in our measured parameters account for the uncertainties due to thermal parameters, modelling parameters, observational systematics, and cosmic variance. Our best-fitting parameters show significant evolution with redshift such that $\lambda _{\rm mfp,H\, \small {I}}$ and $\langle f_{\rm H\, \small {I}} \rangle$ decreases and increases by a factor ∼6 and ∼104, respectively from z ∼ 5 to z ∼ 6. By comparing our $\lambda _{\rm mfp,H\, \small {I}}$, $\langle \Gamma _{\rm H\, \small {I}} \rangle$ and $\langle f_{\rm H\, \small {I}} \rangle$ evolution with that in state-of-the-art Aton radiative transfer simulations and the Thesan and CoDa-III simulations, we find that our best-fitting parameter evolution is consistent with a model in which reionization completes by z ∼ 5.2. Our best-fitting model that matches the $\tau _{\rm eff, H\, \small {I}}$ distribution also reproduces the dark gap length distribution and transmission spike height distribution suggesting robustness and accuracy of our measured parameters.

 

The mean free path of ionizing photons,λ_mfp, is a critical parameter for modeling the intergalactic medium (IGM) both during and after reionization. We present direct measurements ofλ_mfpfrom QSO spectra over the redshift range 5 <z< 6, including the first measurements atz≃ 5.3 and 5.6. Our sample includes data from the XQR-30 VLT large program, as well as new Keck/ESI observations of QSOs nearz∼ 5.5, for which we also acquire new [Cii] 158μm redshifts with ALMA. By measuring the Lyman continuum transmission profile in stacked QSO spectra, we find${\lambda }_{\mathrm{mfp}}={9.33}_{-1.80}^{+2.06}$,${5.40}_{-1.40}^{+1.47}$,${3.31}_{-1.34}^{+2.74}$, and${0.81}_{-0.48}^{+0.73}$pMpc atz= 5.08, 5.31, 5.65, and 5.93, respectively. Our results demonstrate thatλ_mfpincreases steadily and rapidly with time over 5 <z< 6. Notably, we find thatλ_mfpdeviates significantly from predictions based on a fully ionized and relaxed IGM as late asz= 5.3. By comparing our results to model predictions and indirectλ_mfpconstraints based on IGM Lyαopacity, we find that the evolution ofλ_mfpis consistent with scenarios wherein the IGM is still undergoing reionization and/or retains large fluctuations in the ionizing UV background well below redshift 6.

 

We present the Sherwood–Relics simulations, a new suite of large cosmological hydrodynamical simulations aimed at modelling the intergalactic medium (IGM) during and after the cosmic reionization of hydrogen. The suite consists of over 200 simulations that cover a wide range of astrophysical and cosmological parameters. It also includes simulations that use a new lightweight hybrid scheme for treating radiative transfer effects. This scheme follows the spatial variations in the ionizing radiation field, as well as the associated fluctuations in IGM temperature and pressure smoothing. It is computationally much cheaper than full radiation hydrodynamics simulations, and circumvents the difficult task of calibrating a galaxy formation model to observational constraints on cosmic reionization. Using this hybrid technique, we study the spatial fluctuations in IGM properties that are seeded by patchy cosmic reionization. We investigate the relevant physical processes and assess their impact on the z > 4 Lyman-α forest. Our main findings are: (i) consistent with previous studies patchy reionization causes large-scale temperature fluctuations that persist well after the end of reionization, (ii) these increase the Lyman-α forest flux power spectrum on large scales, and (iii) result in a spatially varying pressure smoothing that correlates well with the local reionization redshift. (iv) Structures evaporated or puffed up by photoheating cause notable features in the Lyman-α forest, such as flat-bottom or double-dip absorption profiles.

 

Ly α forest decomposed into Voigt profile components allows us to study clustering properties of the intergalactic-medium and its dependence on various physical quantities. Here, we report the first detections of probability excess of low-z (i.e z < 0.48) Ly α absorber triplets over redshift-space scale of r∥ ≤ 8 pMpc (Mpc in physical units) with maximum amplitude of $8.76^{+1.96}_{-1.65}$ at a longitudinal separation of 1–2 pMpc. We measure non-zero three-point correlation ($\zeta = 4.76^{+1.98}_{-1.67}$) only at this scale with reduced three-point correlation Q = $0.95^{+0.39}_{-0.38}$. The measured ζ shows an increasing trend with increasing minimum H i column density (NH i) threshold while Q does not show any NH i dependence. About 88 per cent of the triplets contributing to ζ (at z ≤ 0.2) have nearby galaxies (whose distribution is known to be complete for ∼0.1L* at z < 0.1 and for ∼L* at z ∼ 0.25 within 20 arcsec to the quasar sightlines) within velocity separation of 500 km s−1 and median impact parameter of 405 pkpc. The measured impact parameters are consistent with majority of the identified triplets not originating from individual galaxies but tracing the underlying galaxy distribution. Frequency of occurrence of Broad-Ly α absorbers (b > 40 km s−1) in triplets (∼85 per cent) is factor ∼3 higher than that found among the full sample (∼32 per cent). Using four different cosmological simulations, we quantify the effect of peculiar velocities and feedback and show that most of the observed trends are broadly reproduced. However, ζ at small scales (r∥ < 1 pMpc) and its b-dependence found in simulations are inconsistent with observations. This could either be related to the failure of these simulations to reproduce the observed b and NH i distributions for NH i > 1014 cm−2 self-consistently or to the wide spread of signal-to-noise ratio in the observed data.

 

We present a new investigation of the intergalactic medium near reionization using dark gaps in the Lyβforest. With its lower optical depth, Lyβoffers a potentially more sensitive probe to any remaining neutral gas compared to the commonly used Lyαline. We identify dark gaps in the Lyβforest using spectra of 42 QSOs atz_em> 5.5, including new data from the XQR-30 VLT Large Programme. Approximately 40% of these QSO spectra exhibit dark gaps longer than 10h⁻¹Mpc atz≃ 5.8. By comparing the results to predictions from simulations, we find that the data are broadly consistent both with models where fluctuations in the Lyαforest are caused solely by ionizing ultraviolet background fluctuations and with models that include large neutral hydrogen patches atz< 6 due to a late end to reionization. Of particular interest is a very long (L= 28h⁻¹Mpc) and dark (τ_eff≳ 6) gap persisting down toz≃ 5.5 in the Lyβforest of thez= 5.85 QSO PSO J025−11. This gap may support late reionization models with a volume-weighted average neutral hydrogen fraction of 〈x_{H I}〉 ≳ 5% byz= 5.6. Finally, we infer constraints on 〈x_{H I}〉 over 5.5 ≲z≲ 6.0 based on the observed Lyβdark gap length distribution and a conservative relationship between gap length and neutral fraction derived from simulations. We find 〈x_{H I}〉 ≤ 0.05, 0.17, and 0.29 atz≃ 5.55, 5.75, and 5.95, respectively. These constraints are consistent with models where reionization ends significantly later thanz= 6.

 

ABSTRACT We compare a sample of five high-resolution, high S/N  Ly α forest spectra of bright 6 < z < ∼6.5 QSOs aimed at spectrally resolving the last remaining transmission spikes at z > 5 with those obtained from mock absorption spectra from the Sherwoodand Sherwood–Relics simulation suites of hydrodynamical simulations of the intergalactic medium (IGM). We use a profile-fitting procedure for the inverted transmitted flux, 1 − F, similar to the widely used Voigt profile fitting of the transmitted flux F at lower redshifts, to characterize the transmission spikes that probe predominately underdense regions of the IGM. We are able to reproduce the width and height distributions of the transmission spikes, both with optically thin simulations of the post-reionization Universe using a homogeneous UV background and full radiative transfer simulations of a late reionization model. We find that the width of the fitted components of the simulated transmission spikes is very sensitive to the instantaneous temperature of the reionized IGM. The internal structures of the spikes are more prominent in low temperature models of the IGM. The width distribution of the observed transmission spikes, which require high spectral resolution (≤ 8  km s−1) to be resolved, is reproduced for optically thin simulations with a temperature at mean density of T0 = (11 000 ± 1600, 10 500 ± 2100, 12 000 ± 2200) K at z = (5.4, 5.6, 5.8). This is weakly dependent on the slope of the temperature-density relation, which is favoured to be moderately steeper than isothermal. In the inhomogeneous, late reionization, full radiative transfer simulations where islands of neutral hydrogen persist to z ∼ 5.3, the width distribution of the observed transmission spikes is consistent with the range of T0 caused by spatial fluctuations in the temperature–density relation.